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SUMMARY 

From September 1976 to September 1979 a meteorogical 
measurment program was carried out in the Sudbury area. The 
objectives of the program were twofold: 

i. to determine the meteorological characteristics most 
affecting the dispersion and dry deposition of locally emitted 
pollutants 
2. to provide a data base for the development of compiiter 

dispersion and deposition models. 
The program consisted of daily releases of pilot balloons and 
minisondes for the measurement of the vertical wind and 
temperature structure of the lower atmosphere. 

After a significant data base was accumulated, it was decided that 
detailed analyses should focus on characterizing the convective 
boundary layer (or mixed layer) in the local area. 

The results indicate that the characteristics of the mixed layer 
play a major role in determining the dispersion and dry deposition of 
local pollutants in the Sudbury area from May through September. 

The average mid-day mixed layer height over the study perioc was 
ICtO m above ground level and the maximum height was 2328 m. The 
mixed layer winds showed an increased frequency from the Jvest, 
Westerly component winds had an average speed of 6 m/s which was 
approximately twice as great as that for the easterly compsnent 
winds. 

The layer above the convective boundary layer capping inversion 
had an average potential temperature gradient of 3.3°C/Km. 
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The average incoming solar flux over the mid-day hours (1100-1500 
hr) averaged for May through September ranged from 600 to 800 
W/m . Those days with known convective activity produced an 
average mid-day solar flux 50 to 250 W/m^ greater than the mean for 
all days. 

The effect of the convective boundary layer on the disperson of 
locally emitted pollutants has been observed to be quite significant. 
The data obtained in this program have been used in the development 
of several computer models to simulate both the dispersion and dry 
deposition of these emissions. 
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CONCLUSIONS 

Data from the Sudbury Environmental Study Meteorological 
Program have been analysed to characterize the behaviour of ihe 
mixed layer in of the Sudbury area. Several conclusions have resul ;ed 
from the program. 

• The mean mid-day mixing height for the months of May through 
September was approximately lO'fO m above ground level. The 
mixed layer height was relatively constant though June, July and 
August and was slightly lower in May and September. 

• The maximum mid-day mixing height measured during the projjram 
was 2328 m above ground level. The maximum for each rrionth 
varied quite significantly. 

• The greatest variability in mid-day mixing heights occurred in the 
months of May and June. 

• Mixed layer winds in the May to September period shox^ed a 
predominant wind speed of i* to 6 m/s. Westerly component winds 
had a frequency roughly four times greater than easterly 
component winds and the wind speeds associated with tl'ie two 
components were noticeably different. For the period from May to 
September, the average mixed layer wind speed for v^esterly 
component winds was 6A m/s while that for easterly component 
winds was 3.8 m/s. The differences in both the frequency and 
speed of the mixed layer winds have major implications :>n local 
air quality and dry deposition. 

• The upper level inversion above the mixed layer capping inversion 
had a mean mid-day potential temperature lapse rate of 3.3°C/Km 
over the study period. 
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• Mean values of total incoming soJar radiation during mid-day hours 
(1100-1500 hr) averaged from May to September were in the range 
of 600 to 800 W/m . Monthly maximum values were 200 to 500 

W/m greater than monthly means. For days with known 

2 
convective activity, the monthly means were 50 to 250 W/m 

greater than the mid-day means obtained using data from all days. 



• The convective boundary layer is a signficant daytime 
meteorological characteristic of the Sudbury area. It occurs on 
relatively sunny days from May through September. Locally 
emitted pollutants released into the convective boundary layer 
behave in a manner known as looping. This is caused by large-scale 
convective turbulence and ultimately results in the rapid mixing of 
the emissions throughout the entire layer, and high ground level 
concentrations near the stack. Because the average mixing height 
is over two and one half times the height of the INCO 381 m stack, 
dispersion of emissions from the tall stack is quite sensitive to the 
existence of the convective boundary layer. 
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1. INTRODUCTION 

From September 1976 to September 1979 the Air Resources Branch 
carried out a specialized meteorological measurement program in 
the Sudbury area. The objectives of the program were twofold, 
namely: 

1. To provide information on the meteorological characteristics 
of the Sudbury area in order to evaluate their impact on the 
dispersion and dry deposition of locally emitted pollutants. 

2. To provide a meteorological data base for the development of 
computer models designed to assess the fate of locally 
emitted pollutants. 

The program consisted of the measurment of wind speed, wind 
direction and temperature in the vertical using a pilot-balloon 
minisonde system. The balloon releases were carried out twice daily 
- once in the early morning and once in the early afternoon -under 
fair weather conditions. The balloons were tracked using dual 
theodolites. Data analysis was carried out in Sudbury using a 
Hewlett-Packard 9825 computer. During the last year of the 
program, total incoming solar radiation measurements were also 
made as part of this program. 

A preliminary assessment of the meteorological data base 
indicated that the program should focus on characterizing the 
meteorology of the Sudbury area during summer months. The 
rationale for this is based on the dispersion of local emissions in the 
convective boundary layer (CBL) which is formed during these 
months (May through September). The convective boundary layer 
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is defined as that portion of the atmosphere affected by heating of 
the earth's surface. This heating results in iarge scale convective 
turbulence which, in turn, causes vigourous vertical mixing of the 
atmosphere. The term 'mixed layer' comes from this characteristic 
of the CBL and is used interchangeably with it. The mixed layer is 
generally limited to the first two kilometers above the earth's 
surface. 

Pollutant plumes emitted into the convective boundary layer 
behave in a characteristic manner described rather graphically as 
'looping'. This behaviour is caused by portions of the plume being 
emitted into the large scale turbulent updrafts and downdrafts of the 
convective atmosphere. Under these conditions, portions of the 
plume can be brought to the surface very close to the stack resulting 
in high ground level concentrations. The existence of the convective 
boundary layer and the looping of plumes has long been noted in the 
Sudbury area. It appeared to be especially important for the INCO 
381 m stack which is the major emitter in the area since in general 
during convective conditions in the summer months the plume from 
this stack was observed to impinge most frequently at ground level. 
Thus, the impact of the emissions from this source on both air quality 
and dry deposition in the Sudbury area is related to the formation of 
the convective boundary layer. It was therefore decided to 
concentrate the data analysis on producing a statistical analysis of 
the convective boundary layer in Sudbury during the summer months. 

This report describes the field measurment program, the data 
analysis methods and the analysed characteristics of the mixed 
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layer. While the statistical results are sometimes referred to as a 
climatology in the body of the report, it should be pointed out that 
the time period and the frequency of measurements were insufficient 
to make this a complete climatology and care should be used in 
extrapolation. 

A detailed description of the diurnal behaviour of the mixed 
layer and its impact on local air quality is published elsewhere (1). 
Two intensive field studies were carried out in the Sudubry area 
precisely for that purpose. The data obtained in those studies 
provided an extremely useful base for guiding the analysis presented 
in this report. 
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2. FIELD MEASUREMENT PROGRAM 



The meteorological field measurement program began in 
September 1976 and continued until September 1979. It consisted of 
two components, namely, the measurement of wind speed, wind 
direction and temperature in the vertical and the measurement of 
surface solar radiation. Both measurements were designed to provide 
a statistical data base from which climatolcgical characteristics of 
the convective boundary layer in Sudbury area could be determined. 

The two measurement programs are described below, 

2.1 Vertical Wind and Temperature Measurements 

A pilot-balloon minisonde system was used for the measurement of 
wind and temperature in the vertical. The instrumentation consisted 
of the following! 

Aero Aqua 30 gm pilot balloons 

Aero Aqua Minisonde Transmitters 

Aero Aqua Model AAl 3001 Minisonde Receiving Station 

\^arren-Knight Pilot Balloon Theodolites 

Sangamo Minisonde Transmitters 

Climatronics Mini-Met Sonde (MMS) System 

Weathermeasure BM70 Micro Aneroid Barometer 

Motorola FM transmitter/receiver and walkie-talkies 

The measurements were done as follows: 

a) A pilot balloon was filled with helium to a calibrated lift of 3 ms"^. 

b) A minisonde was attached to the balloon at the end of a 3 m length 

of string. 
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c) The sonde was adjusted to transmit at a central frequency of '*03 

MHz. 

d) The ground level temperature near the sonde was measured (using a 

digital thermometer) and noted on the chart recorder output 
of the minisonde receiving station. Ground level pressure 
was also measured and noted. 

e) Observations were made on the cloud cover, the weather, 

and the INCO 381 m stack plume behaviour and direction. 

f) Two Warren-Knight Pilot Balloon Theodolites were set up at the 

ends of a 300 to 'fOO m baseline oriented approximately 
perpendicular to the wind direction. 

g) The two theodolites were levelled and sited on each other. The 

angle above or below the horizontal {i.e. the ground plane 
angle) was measured and noted by each theodolite operator. 

h) Both theodolites were sighted on the balloon at its position before 
release. These initial azimuth and elevation angles were 
recorded on cassette tape recorders. 

i) The balloon was released. Every thirty seconds during the ascent, 
the position of the balloon in space was marked by measuring 
the instantaneous azimuth and elevation angles from the 
theodolites and recording them on tape. Also, at 15 second 
intervals, a timing mark was placed on the chart record 
output of the minisonde receiver. From 1976 to 1977, the 
Climatronics Mini-Met Receiver was used with Sangamo 
sondes. This system required an operator to keep the chart 
record ai scale, to mark the 15 second time intervals on the 
chart record and to broadcast the 30 second timing marks to 
the theodolite operators. In 1978, the Aero Aqua system was 
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implemented. The timing requirements were automated with 
this system and no operator was required for the minisonde 
record. The entire release operation changed from a three-to 
a two-man operation, 
j) The balioon was sited every 30 seconds until it was no longer in 

view of both theodolites. The minisonde signal was recorded 
until it was lost. 

Balloon releases were carried out on weekdays only (except during 
intensive field studies) under non-precipitating conditions. Under 
ideal conditions, two releases were carried out daily ~ one in the 
early morning immediately after sunrise to measure nocturnal 
conditions and the second in the early afternoon-1300 to 1*00 hours 
local time to measure daytime conditions. 

The location and details of the it baselines used throughout the 
study are indicated in Figure 1. All were located within 5 Km of the 
INCO 381 m stack. Baselines 3 and tt were most frequently used 
during the program. 

Figures 2 to if show the instrumentation used in the measurement 
program. 

§M Solar Radiation Measurements 

From June, 1978, to September, 1979, measurements were made of 
total incoming solar radiation. A Kipp and Zonen Model CM5 
Soiarimeter with a wavelength range from 300-2500 nm was used. 
The voltage output was recorded directly on a strip chart recorder. 
Before being placed in the field, the soiarimeter was 
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calibrated at the Radiation Laboratory of the Atmospheric 
Environment Service In Downsview, Ontario. 

The solarimeter was mounted on top of the pumphouse of the 
Water Pollution Control Plant on Kelly Lake Road in Sudbury. This 
location is immediately adjacent to baselines ^ and 5 used in the 
Pibal-minisonde program. 

The 1978-79 solar radiation measurements were augmented by 
similar measurements made by the Water Resources Branch in the 
summer of 1977. 
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3. DATA ABSTRACTION AND PRESENTATION 

The data analyses performed to abstract the raw pilot baJloon- 
minisonde data into vertical profiles of wind speed, wind direction 
and temperature are described in this section. The abstraction 
methods for solar radiation data are also described. These analyses 
were the first step in preparing the data for statistical determination 
of convective boundary layer characteristics. 

3.1. Pilot Balloon - Minisonde Data 

Abstraction of the pibal-minisonde data from raw to final form 
consisted of four steps: 

1. Raw data transfer to a Hewlett-Packard 9825 computer 

2. Data analysis using Thyer's method (2) 

3. Re-analysis using data smoothing 

'*. Printing, plotting and storage of profile information 
Each step is described separately. 

3.1.1 Data Transfer 

After each balloon release, the raw data were transferred from 
tape recordings and chart records to summary tables. These raw data 
consisted of the azimuth and elevation angles of both theodolites 
taken every thirty seconds and the ambient air temperature taken 
every 15 seconds. Additional observational data were also included. 
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The information on these summary tables was transferred to a 
Hewlett-Packard 9825 computer. Entry was done manually at the 
keyboard. Data storage was originally done on cassette tape but was 
changed to disc mid way through the program. As the data were 
being transferred, they were carefully screened by the operator to 
eliminate obvious errors, i.e. inaccurate or questionable data were 
either changed to reasonable values (if a straightforward error was 
made) or eliminated completely. 

Table 1 is an example of the tabulated raw data computer output. 
It contains the following: 

a) Date: day/month/year of the release 

b) Time: local time of the release 

c) Baseline: number, bearing (with respect to true north) and length 

(in metres). 
d)Met Data: wind and cloud conditions 

e) Plume data: behaviour and direction of the 381 m stack piume 

f) Surface Pressure: millibars 

g) Relative Humidity: percentage 

h) Column 1: time of the theodolite readings after balloon release 

(minutes) 
i) Columns 2 and 3: the azimuth and elevation angles of theodolite 

number 1 at each 30 second timing mark (degrees) 
j) Columns ^ and 5: the azimuth and elevation angles of theodolite 

number 2 at each 30 second timing mark (degrees) 
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k) column 6: the temperature of the ambient air at each 30 second 

timing mark ( degrees Celsius) 
1) Column 7: the mid-point temperature taken at the 15 second mark 

between theodolite readings (degrees Celsius). 

3.1.2 Thyer's Analysis 

The abstraction of the azimuth and elevation angles to vertical 
profiles of wind speed and direction was carried out using the method 
proposed by Thyer (2). For details of the method, the reader is 
referred to the original publication; however, a brief description is 
given here. 

Thyer's method uses all four theodolite angles to make an optimum 
estimate of the balloon's position at each measurement interval; it 
also indicates the probable magnitude of error in its estimates. 

The geometry of the double theodolite system is shown in Figure 5. 
It can be seen that the balloon's position in space can be defined by 
only 3 variables; however, the two theodolites provide U variables 
(angles) for this purpose. Thyer's method takes advantage of this 
extra piece of information in calculating the balloon's position. 

Referring to Figure 5, one can draw rays from the two theodolites 
towards the balloon according to the measured azimuth and elevation 
angles, but the two rays may or may not intersect. In practice, these 
rays selddm do intersect because of levelling, orienting and 
collimating errors in the theodolites and/or unsynchronized readings 
by the two operators, errors in reading the angular scales, and 
difficulty in pinpointing the balloon position at the timing mark. 
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Even in the absence of these errors, the limitations of accuracy in 
the theodolites can result in non-intersecting rays. For example, if 
the readings are made to the nearest 0.1 degree, then the balloon 
could be anywhere within a cone of semi-vertical angle 0.0^ degree 
whose axis is the ray from the theodolite to the position detennined 
by the measured azimuth and elevation angles. Using double 
theodolites, the balloon could be anywhere within the region of 
intersection of these two cones whose axes are defined by the 
theodolite angles. One can therefore only calculate a most probable 
estimate of the balloon position within the region of intersection of 
the cones. 

Thyer uses a vector method to find the line representing the 
shortest distance of approach between the two rays. He then 
subdivides this line (known as the 'short line') in the ratio of the 
diameters of the cones where it intersects them. This point is taken 
to be the most probable position of the balloon. 

In Figure 5, rays AB and DC represent the cone axes from 
theodolites 1 and 2 respectively. Line BC represents the line of 
shortest approach between the two cones. By defining each line as a 
vector, t.e. AB as r ^ (a^ + bi + ck), DC as r ^ (di + e^ + fk), and BC as r^ 
(gi + P2 + qk) it IS possible to derive the coordinates of the most 
probable position of the balloon as: 

X -r J a+ Pr , /(t J +r 2) g (1) 

y =r J b + ^jlj /(l J +X2) P (2) 

z =X 1 c + j.r| /(r ^ +I2) q (3) 

where 

g = a/A (i,) 

P = ^/a (5) 

q = 7/A (^) 
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and 

a ^ sin Ej cos E2 cos A2 -cos Ej cosA^ sin £2(7) 

^ = cos Ej sin A J sin E2 -sin Ej cos E2 sin A2(8> 
7 = cos Ej cos A I cos E2 sin A2 -cos E, sin A. cos E2 cos A- (9) 
A = (a^ + 0^ ^ y^/^ (10) 

E| and E2 = elevation angles of theodoJites ! and 2 respectively 

A^ and A2 = azimuth angles of theodolites 1 and 2 with respect to line 

AD 

In addition to the coordinates x, y and z for each double theodolite 
reading, the method also provides a measure of the relative error in 
the balloon position: 

r= £ /(r , + I 2^ ^^'^ 0-0^^ ii® 

This, in fact, represents the ratio of the actual uncertainty to the 
theoretical uncertainty (+0.05 degree). 

A computer program was developed in 1977 to analyse the raw data 
according to this method. The program was continually upgraded 
throughout the study and is documented in several internal Air 
Resources Branch reports (3),('t),(5),(6). The program used Thyer's 
method to calculate the coordinates of the balloon (x,y) at each 30 

second interval into the flight. It then calculated wind speed, V , as 

xy 

the vector difference between two consecutive readings divided by 30 
seconds, i.e. 

V=<^\*4/'/30. (12) 

The wind direction was determined as = tan ~ (V /V ) with the 

X y 

following conventions: 
(i)f or V^ > and V > 0, 

wind direction = bearing to true north {N) + 6 + 180 (in degrees). 
(ii) for V^<: Oand V > 

wind direction = N + ^ + 180 
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(iii for V :> and V < 

X y 

wind direction - N + •* 



(iv) for V < and V < 

X y 



wind direction - N + ^ 

The wind speeds and directions calculated in this manner therefore 
represent averages over the height interval between two consecutive 
balloon positions (approximately 90 m). For each calculated wind 
speed and direction, the program also calculated the relative error 
using equation (II). 

For times during balloon ascents when one or both readings from 
one theodolite were missing i.e. only single theodolite data were 
available, the program used a different analysis technique. This 
technique involved estimating the height, z, of the balloon a the time 
of the missing reading and then calculating the other two coodinates, 
(x and y) from the single theodolite data. 

The value of z was calculated as plus 30 seconds times the rise of 
the balloon in the previous interval. Two cases were considered, i.e. 
(i) Theodolite 1 data availble: 
z estimated 

X = z (sin A /tan E.) (13) 

y = z (cos A J /tan E^) ([4) 

(ii) Theodolite 2 data available: 

z estimated 

x=z (sin A^/tan E^) (15) 

y»z (cos A-/tan E2) (16) 

Before the next stage in the data abstraction was initiated, the 
results of Thyer's analyses were manually screened. If inordinately 
large values of r (Eq. 11) were observed or the profiles looked 
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unreasonable, the input data were eliminated (or in some obvious 
cases, changed) and the single method was invoked. 
The temperature data from the minisondes were incorporated in the 
analysis by assigning the proper height z to the appropriate 
temperature measurement. 

Table 2 is an example of the output resulting from the analysis 
using Thyer's method. The following information is included in the 
output: 

(a) Date and Time of Release 

(b) Baseline Number, Bearing and Length 

(c) Time after release (in minutes) 

(d) X coordinate with respect to theodolite 1 (metres) 

(e) y coordinate with respect to theodolite 1 (metres) 

(f) z coordinate with respect to theodolite 1 (metres) 

(g) Height (MSL) = z + 256 metres above mean sea level 

(h) Temperature (°C) at each measured z and at each midpoint 

(i) Rise Rate (m/s) = z/30 seconds 

(j) V (m/s): the interval average wind speed between 

xy 

consecutive readings 
(k) Wind Direction relative to true north 

(1) TAU: the relative error in the measurement. Note that TAU 

is not available for single theodolite calculations. 

In 1979, the decision was made to calculate and present 

potential temperature rather than measured ambient temperature. 

Potential temperature is defined as that temperature obtained by a 

parcel of air at temperature T and pressure P after adiabatic 

expansion or compression to 1000 mb, i.e. 

= TdOQO^^^P (17) 

P 
where R/Cp = 0.286 for dry air. 
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The derivation and method of calculating potential temperature in 
the pibal program is described in Appendix I. It should be noted that 
in order to carry out this calculation, a measurement of the 
atmospheric pressure at ground level at the time of each release was 
required. For historical pibal releases with no pressure data, the 
Atmospheric Environment Service pressure measurements at Sudbury 
Airport were used after being corrected for the difference in 
elevation. 

The decision to change to potential temperature was made because 
it provides a more easily interpreted indicator of atmospheric 
stability than does ambient temperature. This reduced the work 
involved in later analyses quite significantly. 

3.1,3 Data Smoothing 

The wind speed and direction profiles calculated using Thyer's 
method were reanalysed using a mathematical routine to smooth out 
random errors in the theodolite data. The smoothing algorithm was a 
modified version of a smoothing routine proposed by Thyer (7). 

A description of the smoothing routine is given in several 
references (3,5) and the reader is referred to these for details. It was 
developed on the premise that a balloon's ascent is relatively smooth 
in all coordinate axes. This may not always be the case because of 
random turbulent fluctuations in the atmosphere but it is generally a 
reasonable assumption. 

The method takes the calculated x, y and z coordinates of the 
balloon, finds an error bar associated with each coordinate at each 30 
second interval during the flight and then calculates the smoothest 
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path through all of these error bars. The original x, y and z 
coordinates are then replaced by the new smoothed coordinates and 
the wind speed and direction profiles are recalculated. 

3AA Printing, Plotting and Storage 

The smoothed wind speed, wind direction and potential or ambient 
temperature data for each release were printed, plotted and stored 
on the HP 9825 computer. Table 3 shows the final print-out of 
smoothed data (indicated under the title) after the potential 
temperature analysis was initiated. Again it indicates the Date, 
Time, Met Conditions, Plume conditions, Pressure and Humidity. It 
then tabulates the following: 

(a) Time into balloon ascent (min.) 

(b) Height (MSL) in metres 

(c) Ambient Temperature (*^C) 

(d) Potential Temperature (°C) 

(e) Potential temperature lapse rate (degrees C/m of height) 

(f) Rise rate of balloon between readings (m/s) 

(g) Wind speed (m/s) - smoothed 

(h) Wind direction (degrees w/t true north) - smoothed 

The vertical wind and temperature profiles were then plotted. 

Figure 6 shows the plot of the same data tabulated in Table 3. 

Before potential temperature was plotted routinely, both the ambient 

temperature and the dry adiabatic lapse rate (-0.98°C/100m) curves 

were plotted* 

Finally, the smoothed data shown in Table 3 were stored on disc. 

These discs are available both at the Air Resources Branch in Toronto 
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and the Northeastern Regional Offices in Sudbury. The pibal analysis 
program including Thyer's analysis and the smoothing routine is 
attached in Appendix 2. 

The smoothed wind speed, wind direction and temperature profiles 
made up the data base from which the climatotogical analyses for the 
Sudbury area were done. 

3.2 Solar Radiation 

The chart records of solar flux were manually digitized into hourly 
averages. To insure the quality of the data several spot checks were 
carried out against planimeter-determined averages. The comparison 
was excellent. After digitization, the data were stored in the HP 
9825 computer and printed in tabular form. Appendix 3 contains all 
of the hourly average solar flux data obtained in the program. 
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t^. ANALYSIS AND RESULTS 

The wind, temperature and solar radiation data bases were further 
analysed to determine climatologicai characteristics of the 
convective boundary layer. The analyses were aided by a knowledge 
of the characteristics of the convective boundary layer on a diurnal 
basis. These characteristics are presented in a report (1) of several 
intensive field studies carried out in Sudbury. They are described 
here briefly. 

The mixed layer generally grows rapidly on sunny days from 1000 
to 1300 hours local time. By mid-afternoon, the mixed layer height, 
or mixing height, reaches quasi-steady-state conditions and remains 
approximately constant until roughly one hour before sunset. At this 
time the mixed layer rapidly deteriorates as convection is no longer 
maintained. 

At any time, the convective boundary layer is characterized by an 
isothermal potential temperature gradient capped by a strong 
inversion. The wind speed is generally constant throughout the layer 
but increases rapidly above it while wind direction is generally 
variable thoughout the layer and becomes approximately constant 
above it. 

'f.l Mixed Layer Heights 

''.l.l Data Analysis 

The diurnal characteristics of the mixed layer discussed in Kaimal 
et al (7) and Vet et al. (1) were used to produce statistical analysis of 
the the mixed layer height. It was a multi-step process, i.e. 
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(a) . For the period from August, 1976 to .September, 1979, all mid- 

day vertical profiles of wind speed, wind direction arid 
temperature were recorded. Table ^ indicates the exact study 
period. 

(b) A new series of computer plots was generated for the 
Isolated data set. These plots consisted of vertical profiles 
of potential temperature only. For the individual months of 
May, June, July, August and September, the profiles were 
plotted in sequence from 1976 to 1979. The entire set of 
these plots is presented in Apendix i*. 

(c) Each profile in Appendix (iwas analysed for its mixed layer 
depth by looking for the height of the capping inversion. The 
analysis was based on the characteristics of the mixed layer 
described earlier namely, the mixed layer height is delineated 
by a strong capping inversion under which a roughly 
isothermal potential temperature gradient is present. 
When the aforementioned characteristics were not clearly 
discernible, other information was used to determine the 
mixing height. First, wmd speed and wind direction profiles 
were examined for typical convective layer characteristics. 
These characteristics were sufficient for determining mixing 
heights for most of the releases. For those which were still 
unclear, additional information was used, e.g., solar flux, 
cloud cover, morning profiles, and plume behaviour 
observations. 



(d) 
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(e) The mid-day mixed layer height for each release was placed 

in a computer data base. The maximum, minimum, and 
arithmetic mean mixing heights as well as the standard 
deviations, frequency distributions and coefficients of 
variation were calculated for each month by combining 
several years worth of data for the given month. 

iiA.2 Results 

The results of the mixing height analysis are presented in Table 5 
and plotted in Figure 7. They are summarized as follows: 

(a) The June, July and August mean mixing heights were roughly 
equivalent with magnitudes of iQ$it+ ^73, 1072+ 391^ and 1089+ 
311 m respectively. All heights are above ground level. For 
the three months combined, the mean mixing height was 
1082+ if 03 m. 

(b) May and September mean mixing heights were also roughtly 
equivalent with magnitudes of 916 and 908 m repsectively. 
These means were approximately 15% less then those for the 
other months. 

(c) The maximum mixed layer height varied considerably from 
month to month. The greatest value occurred in June at 2330 
m and the least value in September at l'f70 m. 

(d) The minimum mixed layer height also varied considerably 
from month to month with a maximum of 520 m in July and a 
minimum of 2^*5 m in June. The minimum mixing heights 
correspond to days of low thick cloud cover. 

(e) The coefficients of variation for the monthly mixing heights 
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indicate that more variability in the height occurred in May 

and 3une than the other months. 
The frequency distributions of the mixing height data are shown in 
Figures 8 -15. They indicate that the overall data set is normally 
distributed and that roughly 90% of the measured mixing heights 
exceed the 381 m stack height (approximately 680 m VISL). 

^.2 Mixed Layer Winds 

The wind speed and wind direction data were combined with the 
mixing height data to produce an analysis of the mixed layer winds. 

tt.2A Data Analyses 

Each mid-day wind and temperature profile on file in the 
meteorological data base was analysed for its mixed layer average 
wind speed and direction. This was done by converting the stored 
polar coordinates (wind velocity and angular direction) at each 
measurement interval to rectangular coordinates, calculating the 
weighted-average wind components (Equations 14(a) and 14(b)) up to 
the mixing height, then transforming back to polar coordinates. 



H , 

^ =S\i■''^/^^ (I'^a) 

H j 
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where V = east-west wind component 
V = north-south wind component 
Ahj = height intervals between theodolite readings 
H = mixing height. 

This analysis was done for all balloon releases in which the wind 
data extended to at least one half the mixed layer height. Such a 
criterion introduces some bias to the results because of wind 
direction changes with height within the mixed layer. However, it 
was felt that the final statistical analysis of the results was 
relatively insensitive to the inherent errors. 
^.2.2 Results 

Mixed layer wind roses were constructed for individual months 
using the available data from 1976 to 1979 (Figures 16 - 22). They 
illustrate a very strong tendency for southwesterly component winds 
and a general predominance of winds in the south to west quadrant 
with a slightly smaller northerly bias. Mixed layer wind speed 
frequency histograms were constructed (Figures 23 -29). The 

monthly and seasonal plots illustrate a predominant wind speed in the 
't to 6 m/s category. 

The apparent difference in wind characteristics from different 
quadrants led to a more detailed analysis. The data were separated 
into the following categories: 
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(a) all data 

0>)( winds associated with a westerly component, i.e., N to W to S 

(c) winds associated witii an easterly component, i.e., N to E to S 

Table 6 summarizes the mixed layer wind speed data from category 
(a) all winds. The monthly mean wind speeds varied from 5.0 m/s in 
September to 6.7 m/s in June with a May to September average of 5.9 
m/s. 

After separation into directional categories (b) and (c) (Tables 7 
and 8 respectively), the wind speeds show a dramatic change. The 
mean wind speeds from May to September were 6.4 + 3.0 m/sec. for 
the westerly component winds and 3.8 + 2.2 m/s for the easterly 
component winds corresponding to a 1,68 times increase in average 
wind speed. In addition, the frequency of westerly component winds 
was roughly it times that of easterly component winds. 

The dramatic directional dependence of mixed layer winds 
suggested that mixed layer heights may also be directtonally 
dependent. The mixing height data were therefore separated into 
two categories, namely, those having westerly component winds and 
those having easterly component winds, and were further analysed. 
The resultant monthly and seasonal means for the two data sets are 
shown in Figure 30. They indicate that the arithmetic means and 
standard deviations for the westerly and easterly sets of the entire 
data set were 1101 ~ ^0*^ m and 981 - 337 m/s respectively. The data 
were not tested for significant statistical difference. 
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Great care must be taken in assessing the mixed layer wind 
analyses. Unlike mid-day mixing heights (which are relatively 
steady-state and representative of entire afternoons), the wind data 
are highly variable in time. The pilot balloon measurements, being 
essentially instantaneous, do not provide time-averaged data although 
they are analysed here as representative of time averages. In 
addition, the number of data points in some months and for some 
components is very small. Hence, the statistical data given in this 
section should be considered as indicators of tendencies in mixed 
layer wind characteristics rather than as climatological data. 

'^•^ Potential Temperature Lapse Rates 

The change in potential temperature with height (i.e. lapse rate) 
has a characteristic behaviour within and above the mixed layer. 
Within the layer the potential temperature is roughly constant and 
invariant with height corresponding to a zero lapse rate and a neutral 
atmosphere. The capping inversion at the top of layer is usually a 
highly stable narrow layer with a sharp positive potential 
temperature gradient. Above this inversion is a thick stable layer 
with a smaller positive lapse rate. 

Investigators (8) have shown that the lapse rate within the layer 
above the capping inversion is a controlling parameter in the growth 
of the mixed layer because it represents the buoyancy structure of 
the atmospheric layer which convective turbulence must penetrate 
and erode in order for the mixed layer to grow. 
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An analysis of the potential temperature data was therefore 
carried out to determine the characteristic lapse rate of this upper 
layer. 

'^.3.1 Data Analysis 

Each mid-afternoon potential temperature profile in the 
meteorological data base was manually inspected. The layer above 
the CBL capping inversion was noted and the average lapse rate 
within the layer was measured. Only those profiles with a minimum 
vertical extent of 300 to 500 m in the upper layer were analysed. 
Statistical analyses were then carried out on the data. 

'^.3.2 Results 

The maximum, minimum, and mean potential temperature lapse 
rates are shown in Table 9. The monthly means range from 2.8 - 1.2 
C/Km in September to 3.8 - 1.6 °C/Km in May with a seasonal mean 
of 3.3^1.7 °C/Km. 

The frequency distributions of the lapse rate data are shown in 
Figures 31 to 37. The probability plot shown in Figure 38 indicates 
that the overall distribution of the data is normal, 
'f.'* Solar Flux 

The importance of solar radiation to the development and 
maintenance of the convective boundary layer has been noted by 
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many investigators (8,9). The measurement and statistical analysis of 
these data was therefore considered of importance to the objectives 
of the Sudbury Environmental Study. 

itAA Data Analysis 

The hourly-averaged values of incoming solar radiation were 

analysed to produce the following: 

(a) Monthly maximum, minimum and mean values of solar flux 

for each hour of the day during each month. Where more 
than one year's data was available for a given month, all 
years were combined. Special attention is again given to the 
months of May through September because of the convective 
activity during those months. Data from the other months 
has been included for completeness. All days were included 
in this analysis regardless of the weather conditions. 



(b) 



The monthly mean solar flux values for each hour were 
scanned to pick out the peak value. This value was called the 
mid-day monthly mean and was plotted against the mid-day 
absolute maximum value. Both of these values generally 
occurred in the 1200 to IWO hour time period. 



(c) The solar flux data between May and September were then 

analysed in more detail by separating out those days during 
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which convective activity was known to occur. This 
eliminated the effect of overcast and rainy days on the data 
and provided data complementary to the mixing height data. 

i*.k.2 Results 

The monthly maximum, minimum and mean values of solar 
radiation for each hour of all days are given in Appendix 5 and 
Figures 39 to 50. The mid-day mean values plotted against mid-day 
maximum values appear in Figure 51 and the average data for the 
May to September period are given in Table 10, 

The maximum hourly value of solar flux reached throughout the 

2 

entire study period was 1123 W/m . The monthly mean values at mid- 

day, however, were 200 to >00 W/m less than the respective monthly 
maximum values. Mean mid-day values ranged between 600 to 800 
W/m typically. The diurnal cycle of the solar radiation was 
sinusoidal as expected. 

Analyses of the convective solar flux data are presented in monthly 
averages for each hour of the day in Appendix 5. In Figures 52 to 5S 
both the convective data and the all-weather data are plotted 
together. Figure 59 is a plot of the convective versus all-weather 
mid-day mean solar fluxes for each month. The average convective 
data for May through September are given in Table 11, 

The data indicate that, for the period studied, the mean mid-day 

2 
solar flux for the convective days ranged from 650 W/m to 850 

2 2 

W/m . These values are 50 to 250 W/m greater than the mid-day 

means for all-weather conditions. 
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5. DISCUSSION 

The results of the mixed layer height, mixed layer wind, lapse rate 
and solar flux analyses must be considered with some care. Because 
of the short term nature of the program, the number of data points in 
each analysis is quite limited. The results reflect only those years 
during which the measurements were taken and extrapolation beyond 
those years should be done with caution. 

The mixed layer height and mixed layer wind analyses are very 
important in assessing the fate of pollutants emitted in the Sudbury 
area. The consistent mid-day average mixing heights of roughly 1000 
m above ground level for the months of May through September 
indicate that the effectiveness of a tali stack for eliminating ground 
level concentrations is minimal during convective periods. This is not 
the case however, during the night and during winter months when 
nocturnal inversions {shown in the early morning minisonde data) are 
generally below plume level and tend to trap the plume aloft- 
preventing it from impinging at ground level. Therefore, under these 
conditions, dry deposition is expected to be minimal. 

The strong bias for westerly component mixed layer winds suggests 
that the dry deposition field for locally emitted pollutants should also 
be biased. Based on the data analysed here, increased dry deposition 
would be expected in the northeasterly, easterly and southerly 
directions due to increased wind direction frequencies alone. 

The indications (based on a small data base) of stronger wind 
speeds with westerly component winds also has a major impact on 
local air quality and dry deposition. It suggests that, on average, 
ground level concentrations and dry deposition, (from meteorological 
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consi derations alone), should be lower near the stack for westerly 
component winds than for easterly component winds. 

Finally, the convective boundary layer is a very significant 
meteorological characteristic of the Sudbury area. Its existence and 
structure have been well-established and a significant 

meteorological data base is available for assessing its impact on the 
fate of local emissions. 
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Mixed Layer Wind Speed Summary 
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TABLE 9 

Potential Temperature 

Lapse Rate Summary 

(°C/Km) 
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TABLE 10 

Ministry of Che Environment 
Sudbury Environmental Study 

Solar Flux Analysis: 
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Ministry of the Environment 
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Solar Flux Analysis: 
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Fiqure 1: Theodolite Baseline Locations 
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Figure 2: Meteorological Program Instrumentation. The pilot balloon, 

minisonde and theodolite are visible in the foreground, the INCO 
superstack in the background. 




Figure 3: Meteorological Program Instrumentation. A Warren-Knight Theodolite 
being used to follow the pibal-minisonde. 
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Fiqure 4: Meteorological Program Instrumentation -- The Kipp and Zonen 
Solarimeter 



\ 



-REGION OF INTERSECTION OF CONES 




I 

i 



HUNWRIES BF COMES ' 
■*fl kHO VECTOH 



Figure 5: Theodolite - Balloon Geometry for 
Thyer's Analysis (From Thyer(2)) 
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Figure 6: Vertical Profile of Wind Speed, Wind Direction and Potential 
Temperature. 
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Figure 7: Mid -day Mixing Height Summary 
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Figure 8: Mixing Height Freauency Distribution - May 
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Figure 9: Mixing Height Frequency Distribution - June 
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Figure 10: Mixing Height Frequency Distribution - July 
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Figure 11: Mixing Height Frequency Distribution - August 
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Figure 12: Mixing Height Frequency Distribution - September 
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Figure 13: Mixing Height Frequency Distribution - June to August 
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Figure 14: Mixing Height Frequency Distribution - May to September 
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Figure 15: Mixing Height Probability Plot - May to September 
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FIG, 16 Mixed Layer Wind Rose - May 
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FIG, 17 Mixed Layer Wind Rose - June 
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FIG, 18 Mixed Layer Wind Rose - July 
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FIG. 19 Mixed Layer Wind Rose - August 
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Fig: 20 Mixed Layer Wind Rose - September 
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FIG, 21 Mixed Layer Wind Rose June - Auqust 
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Fife: 22 Mixed Layer Wind Rose May -- September 
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Figure 31: Potential Temperature Lapse Rate Distribution - May 
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Figure 32: Potential Temperature Lapse Rate Distribution - June 
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Figure 33: Potential Temperature Lapse Rate Distribution - July 
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Figure 34: Potential Temperature Lapse Rate Distribution - August 
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Figure 35: Potential Temperature Lapse Rate Distribution - September 
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Figure 36: Potential Temperature Lapse Rate Distribution - June to August 
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Figure 37: Potential Temperature Lapse Rate Distribution - May to September 
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Figure 38: Potential Temperature Lapse Rate Probability Plot - May to September 
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Figure 51 
Variation of Solar Radiation Mid-day Maxima and Means 
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COMPARISON OF MEAN SOLAR RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 
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COMPARISON OF MEAN SOLAR RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 
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COMPARISON OF MEAN SOLAR RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 
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COMPARISON OF MEAN SOLAR RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 
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COMPARISON OF MEAN SOLE RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 



r 2013 



r\ 


[ amB 


ni 




<r 




51 




\ 


B0EI 


X 




\J 




X 


E00 


n 




J 




u. 






H00 


in 




n 




_i 




n 


ZBH 


Ul 









1 1 1 1 1 T 1 1 1 1 1 1 T — T — n 1 1 1 1 1 1 r 

UUNE-RUBUST 

+ «=CDNVE:CT I VE C>HY5 

n!=RL!_ DHYS 



' iE. ' . ID ' a r^ 



3 H S B 




H^3 2H 



i 
f 



HOUR 

FIGURE 57 



COMPARISON OF MEAN SOLAR RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 
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COMPARISON OF MEAN SOLAR RADIATION WITH MEAN RADIATION IN CONVECTIVE PERIODS 
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Figure 59: Monthly Mid-day Solar Flux Data 
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APPENDIX 1 



Calculation of Potential Temperature 



Definition: 9=T/1000 \ R/c (A-1) 



T / 100 



From the minisonde flights, T is known as a function of height but P is 
not. Since it is required for the calculation of O , it must be found 
functionally. This can be done by using the hydrostatic equation and the 
perfect gas law, i.e., 



=; - P9 where p =: p/rt 

dZ ' 



dP _ Pg (A3) 

dZ ~ RT 



dF_ _g dZ (AU) 

p ~ Rt 



If one considers the format of the minisonde data, i.e.,T:/^z) measured at 
15 second intervals, then the atmosphere can be represented as a series of 
layers. Each layer has a depth determined by the rise rate of the balloon 
and is considered to have a linear temperature change with height, see 
Figure Al P;^^ ^ T, ^ , 

Layer 

P: 



\- 



I 
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In layer i , Tz = Ti + 7^ where 7 = Ti + , -Tj 



Equation Cf) can then be written: 



'"{'^y-^l^ 



z 
Z: 



H7 



(H 



= -9A 



^Zi in (Ti+A 

i + ,-Ti) \Ti / 



. . -gAZi/R(Ti-H-Ti) 



r) 



[r] 



(A5) 



(A6) 



(A7) 



(A8) 



, 9AZ/R(Ti-Ti+i) 
P: + , = Pi|Ti+i I (A9) 



Using equation {A-9), the pressure at each temperature measurement 



-102- 

interval can be determined. This value of P can then be substituted back 
into equation (J) to determine the potential temperature at each height. 
All that is required for this calculation is a knowledge of the ground level 
pressure. 

The following procedure changes the ambient temperature profile to 
potential temperature within the analysis program. Figure A2 illustrates 
the situation. 



Surface 
\ 



-» ■ ■ 



AZ, 



AZ, 




T„ 



K^C) 



Fig. A-2 



1. Store temperature and height data for each 15 second interval of the 
balloon ascent. 

2. Calculate the surface level potential temperature as: 



♦■*> 



e,= T /^\R/Cp 



(A-10) 
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3. Calculate the pressure, Pjt^ at each temperature point in the profile, Tj+^ 
using Equation (9), i.e.; 



•"i+i =''i 



Ti^gAZ/R(Ti-Ti + ,) 



] 



where i = 0, N and = surface 



(A-11) 



^. Calculate the potential temperature at each point i using Equation (A- 



11). 



^i^.-"^^. 



lOOOl R/C, 
-Pi*1 J 



Note: Where ground level pressures at the time of balloon release were 
not available, they were estimated using Sudbury Airport pressure data 
accounting for the difference in elevation. 
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Appendix 2 

5ES Pilot Balloon - Minisonde Analysis Program 

(for Hewlett-Packard 9825) 
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2: ^iTi U^^f 31 ,a5[ 31 ,: !?[31 ,1[ 34) ,M34| ,'>r341 ,-M34| 

3: iiT\ V$( 3] ,05[ 31 ,T5[121 ,'{$[ &1 

4: "3t^^r':a3o "t3no only? 3nter-39 for bassi ino " ;«?3it 2030 

5: ?nt "T^SPLT^TF: f? " , cO , "GR^a-n 9L'>^]^: MGL^?" , rl, "^'E^^ T\ r,»? " , F ?■[ I ,30| 

6: 3nt "PLaME O'VP'X?" , 3$[ 1 ,33) , "DAY?" , ro4 , ".'■;OMTH?" , r5 5, " YF^''?" , r6G 

7: 3nt "TI.^^E 0? RJ J?" ,r67, "outout?", 1? 

U: 3nt "PPP''-"JR;^?" ,C> "PFL 'TJ -inT FY?" , P 

9: if c 0=1; 335. yH.r37r 68. 23 + r33; 259. 15*r30; jmo 5 

10: it c 0=2; 3;J0. 4*r37r 166. 23 + r33 ; 259. 15>r33; jmo 4 

11: i£ c 0=3; 343. 95*r37; 237. 3*r33; 256+r30; jmo 3 

12: if c 0=4; 325. 'J4*r37; 154. 5*r33; 256>r30; jmo 2 

13: ort "base" ,rO;ent "Base length ?", r37 , "Ffear ing ?", r33 ; 256+r30 

14: if r 0=-99; sfg 8 

15: fut 0,50x, "ONTAPIO MINISTRY :iF Y^\^ ll'^lVO'^'AZ-^V,/ 

16: fTit l,54x, "SuSbury Ehvir omient al Stuiy",2/ 

17: e-nt 9,10x, f2.'0, "/", f2.0, "/", f4.0,15x, f4. 0,4/ 

18: fut 2,3x, "Ti^?",5x, "Azinuth 1" ,3 x, "Els va tion 1 " ,4 x, "Azi nuth 2" 

19: f-nt 3,51x, "F,l3vitioT 2 " ,6 x, "TTJTn. " ,,■> x, ""15-?oint Tgtio ." 

20: L- Tt 4,3x, " (nil )", 4 x, " (a?7r3e=3) " ,4x, " (i^iC333) " ,5y, " (i?^r33.-;) " 

21: fnt 5,52x, " (l3-3C9a9} " ,4x, " (JJ39r22o C) " , S x, " {^ ^gr ^^3 C)" 

22: fnt 5 , 5 x, "PPFS^O^^" , 3 x, f5 . 1 ,10 x, "pt^L . -lUMIlI ry' , 3 x, f4 . , "'^ ", 3/ 

23: Eit 7,5x, ";tet. Oita", c;35 ,2/ ,5 x, " Plune mti",G3G,2/ 

24: f-nt S,44x, "PrsiBlin? ", f 2 ."0 ,3 x, " I'oar i ni ", f7 . 2 ,5 x, "L3ri^tri ",f7.2,3/ 

25: is? ".53t orinter to .5, nresr, ClM''";3t? 

26: .vti 6,27,B4;v/to 6 ,2 7 ,70 , int {1055/64} ,L Tt (10 35) 

27: vt;^ 3 ,2 7,/6, int (95 0/6 4) ,Lnt (950) ;.Jf t 3;wj:t C. 1 

23: vrt b.9,rS4 ,r.i5,r35,r')7; wto C,27,10;wtb C,27,10;wtb .'J,27,10 

29: LE l:L-) 6=0; Afct b. ?i , rO , r33 , ri? 

30: vrc 6. 7 , ^^^,3 ^'; vr t j . 5 , 'J, "?. 

31: ;t:t G.2;.vtT 6 ,2 7 ,1 l; .;r t 6. 3; vet G.4;vtb G,27,lJ;wrt *;; . 5 

i^.: .vto u,27,l.i:-TC 1^1 to )3;'/tb t-,55;'i-xt Tj.vtb f,13;wt> .:i,10;,-?t- 5,1 

n 

33: eit G,lx, l3. 1,:. x, 1:7. 2 , 3 x, f/. 2 ,7 x, -'7. 1 ,'S x, f 7 . 2 , 3:: , ''S . 1,1:,' x, fS. 1 

34: fT.t 2,lx, r3. 1,'i V, f7.-2,-5x, c7. 2 ,7x, t7.2 ,5x, f7. 2 

35: cnt 3,lx, f3. 1, 3 X, t7. : ,5 X, t;7. > ,7 X, f7. ^ , '. X, f 7 . 7 ,3x j*""!. 1 

36: cnt ^ ,1 x, t :> . 1 , j.i k, C3. 1 ,10 x, f i. 1 

37: cnt J ,1 x, f 3 . 1 , 1 J x, t5, i 

33: "3qo" :J3C 1 = 1 to 33 

3r: if T41;if .' 1 i *^;; 2*O*r40; j-i^ 2 

40: ^Tt "^::i.^ut"l 1? ' , d ) , " '"L ?v L ? " , rll , " '^•^i vjf i 2 ? " , r42 , " iT n' 2?",r;3 

41: :>,it "'i-1-^Dint 'rBm .■'" , rn , ""03. 'Joint 'to ■!'; .?" , ri '1 

42: c 4u*A [I ];: ll*MT];r42+:(Il r r ^2*n[l ];r 44 *■■ [I ,11 ;r --J '.-.-:■; f I ,21 

43: " Jki ": (t-l) /2*r23r ic i[I]=2":J and F [I ,21 = :;■;!; 0* 3 [ I ] ; 1 1-. "jr>n:i" 

i4: if ci.3a=0;it afi]=2"9 an^ f [i ,2| H 99: -. f :j 3;T-l*r5a 

4 5: it fljo anJ [ = l;.7rt j , ri 3 ,i; [ I ,21 ; q to "kne" 

46: i£ eiit;;.vrt 0. ^ ,r23 /"^ [ ^ ,21 ,'^. f T ,1] ; ito "kno" 

^7: it I=l;./t:t b. 3 ,r > 3 ,\ [ [ ] , -> [I ] »:[ H ,T [ M ,^ [ [ ,2} ; ] :i-i 1 

4rt: i E '■ [T ,11 iiOO; I'Tio 2 

49: vet o.^,r?3,\[ri ,nfT] ,:[n ,Tm; ini 1 

50: vrt o. 6 , r2 3 ,\ t ^ ] , B [II ,2 [ I ] ,1 [ I ] ,^ [ [ ,2| , U f ,11 

51: " Kn? " :ri3xt I 

52: ''.10.13" :I-l*r J3* T;i I: f I :) 3 ; r TO *r3 5; I *r 50; r 3 5-^T 

53: "six" :"733"*a$;i f r35n;r35*T 



-ice- 



's 4 



■J: 'nt ),^x, "-nnn- " , i:: . J ,5 x, " -^-ir ini " , f 7 . 2 , 5 x, "r,-nit i ",^7 :>,3/ 

"j: it ll ] 2; }r3 "^^l^" 

53: /t:> '',,12 

57: vcc ir.vrt (;. 1; vr t G. '^ , rS 4 , rj 5 , r3 .3 , r3 7 

5J: fit l,lx, "'rii,",/x, ■' i'MDx, "r,.]x , " i)^!^ 1 1 " , 3 x, •■-^i'vit (^-^r, ) " , ^3 x, "I-ti 

S'): Cnt 2,73x, "ri^? ^':i t V , ) x, " Vx /ML x, ",;i n 1 Oi c 3cti ^.V , j x, " lau " 

aO:^ tut 3,lx, " (niT)",3x, •' (natr-T) " ,3x, "(.n9tr33)" ,3 x, " C^^tc-^)" ,3x, "{nst 

TGO) 

61: fnt 4,^.3x, ■■ (j-jr333 C )" , 4 x, " (notr 35/ sec)" , 1 x, " (Ti3tc 33/ sec )" 

6 2: frit j,101 X, " (.',3-3 . :3i. ttrue north)" 

63: ./tD o,13;.JCt f, G , rj , r:3 3 , r3 7 ; Wtb 0,13;wrt 6.1;wto 6,27,in-^rL b.2;wr . 

t b . 3 

^34: vt^ 0,27,10; vrt 6,4rwto 6,27,10;-vrt o.Srwtb 6,27,10 

65: r:,c m = i tc L3J;vto L,S5;i2Xt 'Tjwt^ 6,13;i/tn y,10 

j5: :nt 2 , 1 x, f., . 1 , ? x, f J . 2 ,2 x, f) . 2 , 2 x, H. 2 ,2 x, f j. 2 

''.7: u' Tt 3,3 2x, f J. 2 ,3 X, E5.-2 

■:o: L-it 'Mx, f J. l,2x, E3. >,2x, ^1.2,2x, e).2,2 X, f).2,23x, fl. l,12x, n. i 

C): tnt o,32x, LJ. j ,3x, es. 2 ,7 X, f3. 1 

7U:^cnt 7,iK, o.l,2x, n.2,2x, n.2,2x, n. ^2 x, f J . 2 ,7 >, t . . 1 , 11 x, H . 1 ,] 2 x, 

71: '-: Tt '3,11 ix, 1:3. i , ; X, Fj. 1 



1-' 1 L 



3,lx, rj. L,M;<, t). 2,3 X, 1). 2 ,/ X, f3. 1 

r n 1 ■ \ '^ 



73: L ii -^ r 1! [■2 " 1; ] t^ 3 

71: J*7.;')>rl.i:!rt ,r23 , ^, L fr ] ) , U 1 ,2 1 ; 2 -f 

7 _^ : 1 1 :> " t ? T I -: L " 

7'.:.: 'J3Lr":L::ir'^ = it3T 

/7: i. L •' !T 1 =2" "I; -) t -'tc^lcl" 

7 3: LL MM"^'^''-' inJ 'M' 1 ^'''1'); !t T "5>iv" 

''■:: 7 3 ^ "3in -il -> " 

l:: tto "3^i^" 

2 1: ■' >;i /"jjT- (■.[.-'!), in ( ^ f^ 1 ) * r 3 ; - ^t (Mr] )CD3 ( ^ f ^ ] ) -r '■ ;i i o f ■ [.-1 j.^-, 

.,.': 3 1 MM^ 1 ) -r.; ji- (T f'' ' )- in {c^ [ ^ 1 ) *r. ; r t M"^ f" M'-oi {Cfi'i ) ^r? 

'3: - 3r/-rl r ;*r '; - 3r i-r". ri-rl 1 ; r '^r',- r3 r7 -r' I 

■■-..i: ■'('-' ■^>rn~^frLl"!) *rL^.; r '•/rL2-rL3;r l../rl2*L-i i;r ]Vri2-rL- 

■..5: : i 3 r > f r L ; r L U r 1 3 r I L > r L 3 : r 3 7 r i n (c L ) * rl 7 ; r 2 7 err. (c 1 ) > rl : 

-.': (c I., (r MrL 3- r^ rL 3( )■ ri /; r 7rl 3- r.) r L 1 ) ) / rL3 *r.3 1 

-.7: fr r;(c jrl3- r.^rLli i-riyr c ^ ri 3- r3 rL i ) ) ^ rl -, -r3 ^: (r T r] n-rL /r' 1 ) /ri , ^r 



J J; 



L ^'iri h(c 2or3 I/: r ■"'■: fr2 3i ) r 12-;; r I'^^rif { r 
r 2 'r .f (£ ?,3ri ;/( r 2^. fr2 3) ) r ]■"-;: 
: 2irin(^ {■ ] ); ji (^ f 



r2 •;/ ( r 2';f L'2_3) ) r I ; v/ 



J J : 
> i : 



r:.' 1 



1 ) * : [ '^ 1 ; t: : .^- i r ('- f ^ 1 ) : :> -, f t | ;- 1 ) » r [ ^ 1 
'L: CM-,: (.". r-^i ]^;, (M'M )M[^l;r :74ri>cr;- (2 ^^ J ) : ^ M M mV ' ^^J 
-^: c i I .-;in C^. |2 ] )* •■ I M;: J / tiM (rL) h r3 3r,in (-1 r M )► P l"i 

[ M-:) *^^:[ ^l>.r n ; (-,[-]-/) *2fM^] -r f^ ];(': p ]-.)*-, ^M - ' f - 1 
I ■■ 1 - ' ) * ■'.+ M M - I ^ ] ; (-1 f M - ■■ ) *2f ^ [ ^ 1 » M" ! ; (.[.]-;,* . + , r -^ 1 
Jo: 'zKi^":ii 2.j;;Ll ' ^ i -> ' ( ■: , r ! 7 , i , r^ i , ;■, r 2 .) , I , : , r , r; - ) 
I'j: *lJ ?; ;*rM; -^ i [, I; n> > r^j. -,,,. f . I . (v +f,,^^ / 3 * : ; ^> r ■ - * r :{, 
:U: r •.o/(r ..Ui:2..) tin ( .■ i) +i; (^-i) /2^r2 ) 
i'3: L L L i -I 2; ]ln " t :;xt r"' 

.■■': ■: It 1,2-, ':3. J ,2 x, [ ) . 2 , > x, f ^ : , 2 x, t") . 2 , 2 x, n . 2 , 7 x, ;-3.1 
10,'): i F -'^rlj-ito -'j-ii-i. " 

17 2: ^ c t t . 2 , r2 3 , < , f , :i, li +: 3 3 ; j ^^ 2 

1J3: vrt 6. I,r23, w^ ,^, 2+:3J,-[? ,2] 

12 1: Tt:> "T^xtF" 

1'35: ■Jv-i":ij: f f"^ , 1] .: D3; j -10 2 

IOj: /: t f . 3 ,./ ,./*-r3 ); j ti 2 
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107 
lOJ 
lUJ 
110 

111 

112 

113 

111 

115 

116 

117 

113 

119 

120 

121 

122 

123 

124 

125 

126 

12 7 

12>3 

12 9 

130 

131 

132 

133 

134 

135 

135; 

137 

133 

139 

110 

141 

142: 

2] 

143; 
144; 
14 5 
lib 
14 7; 

14 3; 
U'J- 
1-jJ 
151 
152: 
153: 
154: 
155: 
153: 
157: 

15 J: 
159: 
lou: 
lU: 
i'o2: 



Ent 
cor 

rfC t 

wr t 



Z$=U$;gta "five" 



vet 6.6,'J,^ + r3),^, [ '■ ,1] 

iC nfr:- ,2] 1*^9; j TiO 2 

vrt 6. 4 ,r23 ,< ,f , Z, Z+i:3 1 ,1 ,3; jmo 2 

.vet G. 7,r2 3,< ,<: ,Z, Z+i:31,'^ [? ,21 ,1 ,3 

vto 6,27,10; rfct C.5,T,l 

":i2xtF" :i£ -^ = 1; X+r47? 7*r43;Z*r A? 

cnt 7,lx, f3.1,24x, f9.2 ,2 X, 0.2,7 X, f3. 1 

i£ ei.j e=0; JTi-^ C 

"tcalcl": 

f-nt 7,lx,f5.'l,2 4 X, rj.2 ,2 X, f?.2 ,7 X, f5. 1 

6 ,31x, fT. ^ ,2 X, f9. 2 ,7 X, f5. 1 

H=P to r 50; CT-l) /2*r2 3 

6. 6, Z + 40. 5,^+40. 5+c3 3,^ ['I ,n 

6. 7,r2 3,Z + 81,2+8H-c3 3 ,S[H ,2) ;Z +81*Z 
113 xt H 

if flg2;gto "-jone" 
"iataolc" :3nt "Data ok?",Z$;if 
cfi 0;q3b "fudge" 
gto "six" 

"fLV3":wtb 6,12;gto "gone" 
"single" :9999999*r26; if F =2; 2. 7-*i; jmo 2 
R[F-1] -»H 

z + 30H+z;if c[Fi ttgggjj-np 2 

Zsin (A[F] )/tan(3 [F ] )+X ; Zcos (A [F ) }/ tan (8 [F ] ) + Y; jmp 2 

Zsin {C[F] )/tan(D[F] )*^;r 37+Zcos (C[F ] ) / tan (D [F ) ) + Y 

X-.005X+K[71; X+. 005X *L [F 1 ; Y-. OOSY-'l f F ] 

y + .005Y+NfF];Z-.005Z+OrFlr Z+. 00 5Z*P[F];cet 

"Euig3":2nt "Time of reading? {-1 to ouit )" , r31 ; if r31=-l;ret 

ent "ThBolDlites or Ifetno .? (1 or2)",r32;if r32 = 2;gto +6 

eat "'ri30lolite 1 or 2?" ,r32;if r32=2;gto +4 

ent "\ziTiuth 1?" ,^ {r 31*2 + 1] , "nlevati :)n 1 ? " , P [r 3 1* 2 + 1] 

if A [r 31*2 + 1] =2"9; r 31*2*r35 

7 to "f'ji33" 

3nt "^^i^lut^ 2?" ,C [ c 31*2 + 1] , "PlsvatioT 2 ?" ,7 [r 31*2 + 1] ; ^to "fuJi?" 

3nt "Mij-o^iat T?-np .?" ,'^[r 31*2+1,11 ,"0:)S. "oint Tenn .?■■ ,i^ [r 31*2 + 1 , 

"313 " :o l-p2+n2r .:> 3-o4 *r>\; o 5-i:^*p5; r> ?/30-*pl4;n 4/30+ra5rT G/3J*n7 

";vnJ " : •/{ p2"2 + p4''2i /3J *?]; 3tn ( p2/r-4) +Pl3; o l^'p2; o 3+rl; ^ 5+oo 

■jtD "/iaV*;if pl4>U and nl3>J or oLKJ ,anl ol3>l;3to "-vinJl" 

" vinJL " :o lU + oi'6 + 13 J+d>; gto "thre3 " 

'" vinj" ::) 10 + 013*03 

"tir33 " :if o V>36 0; o 9-36'J*p; 

if :» y<J; 3 ^+36 0+pJ 



: t: 3. 



jfg /;itD 



.iioth" 



" 1^13 " : 4 6*r2 2 

if A[l]=2"-* 

if jc 1=1 to r33 

Tax (KfTl ,L[T] )»i: ); nin ('-[Tl ,L\ 

nix (MfT] ,1 [II ) +r3; iiin (-'fl] ,Tf 

■l^'< (^(I] ,P[T 1 l^rl; nin O [J ] ,Pf 

1 3 ;< t T 

c :>c 1=1 t^ 3 

c 22+T*r 33;rr 33*:{ 1] ; l*'^ 

"^^^^"rif: P>=:.13?7to "Dut" 



I !)^ 
M ) 
tl )' 



h\i 1 

c'fT ] 



: ) * U T 1 
: 1 * If T 1 
: l + '^ri 1 



:[^1 



*-c 3:*;'^ +1+1? 
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13 3: 


1j4: 


i.:o: 


LjJ: 


La 7: 


10 3: 


1 J D: 


17.J: 


171: 


172: 


173: 


171: 


1 75: 


17j: 


177: 


173: 


17 0: 


13 U: 


lU: 


1J2: 


i:3: 


13 4: 


i:.!.: 


1j5: 


lil: 


13 3: 


IBS: 


iJJ: 


I'Jl: 


l';2: 


1^3: 


1 >1: 


1^3: 


I-Ij: 


1M7: 


]'j: 


}-j^J: 


2 13; 


kLU: 


J I'll 


233: 


;: ■ -1 : 


r;j5: 


LiiG: 


z..i7: 


: ,1 3 : 


tic:- 


.:Ui: 


^^7,1 


3 11: 


212: 


^^. 1 


^13: 


Jl^.: 


:;j5: 



(■■aM-n/( (-'^33) *3J) *L'[-M; (r. [-•'!-.;) /({ ,^-.-.33) *3J) *vrFl 
If: '^+1>.:33; : f ( J r ^ ] -ny [^ l ) * 3 i/2 ->•: f ^ 1 n ho "out." 

{•{[• +11 - :) /( { -^ fx-r35) *33> *UC^+11 ; ('. f ^ + 1| - 7) / { ( ^ +l-c3'>) *3'l) -V f^ ^ 1] 
in vMM >J [^' +11 ; ^ [M *■: [ ^1 Mto "oack" 

■^+l*'^;ni:i (.1 r^l ,J ft^-n )>'Jf^l;Tl^X (V [^1 ,1 \^^ -1\ ) * V \? ] 
if '•'+1<=: 33; it:* "tT3re " 

?Hj( )[=■) +/f M ) +3[ '■-H ; J + 33( J f^l+\/ r^M ) •*"[^]; ito "out" 
■' t -1 3 c 3 " : ( :: [ ^ + ] 1 - J) / ( ( ^ + 1- r3 3 ) * 3 1 ) + 1^'=' + 11 
('3 r •' +11 -!)/(( ■ +1-: 3 3) * 31) *Vf!='+ll 
it" / f^^+l! <-,M-' 1 ; It") "T032" 

if J f?] =T r-"-!] ; K[1-_11 *-f 7-1) . c_l+F.r,tO "013k" 

i + 1 [M *3J*:[^-11 ; K[^] *:r ^1; 7to "back" 
"-J33?.":if y [?+ll >-/ [^l ; ito "oo-, 3" 
in / [?-n =/ [^]; b ['^-11 +:[ ^-11 ;?-l+P;ato "oi-k" 

?+■■/ [^M *3J*:r ^ -11 ; r,[r] +3t^l Mto "b^ck" 
"OT' 3" :?-l + '?; 
" 1 n^" :'^+l* "^ 

( M ■ !-?)/( (^-:33) *33) *nfPh (r.t^i- ;)/( {• -:.33) *30) + V fp ] 

if / [l^I <=J [■='-1] inl t.i ["^l >=7 ['^-n ; no "loo" 

? - 2* c 3 3 

"Ijc " :cr>r X=Z to c 20 

v.' + { J [^-1] +; [^-11 ) *15*(X- r33» *:;[ <\ ;nc!Xt X 

K-l*Fnto "•:iJZ<" 

"lDo":Tiin (-J [?] ,J [?-l] }*u [P1;ti3x {\/[?] ,7 [^-11 )*V[F] 

if ^+1>!:33:""' +l+F;F-l+r23;ito "for" 

7 to "agaii" 

"o'jt":if I<2;ar3 C*n,;ir3 l+.<;ara "*r.rJ!T'n 2 

i I; I<3;a£:3 C+;"far3 0^;ar3 p* L 

13 xt I 

it" : L J 3; 3 t o " ; t " 

"Bnotti" :it flj 7; 2*.C; 0-3[ 11 ; j-ro 2 

i£ Ll7 C;r 35+l+< 

it ':Li7 or £173; for J =K to 3 1 ; C[ J -1] +81. 5+2 [ J 1 ; n3xt J 

if >"«"/33" JlTtO "3 0" 

cut 'j,5ijx, ".'rrrA no ^'T.nsr'^Y ^^ ri3 3 i/t ?") n-i 1 ["i, / 

-.It 1,54 X, "^al :iat: V ^n/ir: oonont^l ''t j 1 y", / ,31 x, " ^n Dot lo 1 -^^ti",/ 

I ,it 2, jx, "Ti ij ^',3 X, "M^iTit (■'^'3 )", 2 X, "T^no", 1 x, "Pot^nti iL To 10 ■• 

r nt i, 3X, "p-"'^:i'j-r",3 X, t3. 1,1') X, "^'^L. . 3i"T:tr ry' , 3 x, fi. j, "■?: ",2/ 

lint %3x, " (nin )", 1 x, " (n^tr-^) " ,2x, " (i 3i ?)" , 1 x, " (i?i r)" 

jTt ->,lSx, "(Km ^/-i) ■' ,^x, "^/r:)" ,5x, " (n.-tc?3/ sec}" 

^ nt ),i6x, "(Kj C3l. tcuo north)",/ 

:it ;-;,'::x, "■ :^r ^^ r>: ' , c3 j ,1/ ,3 x, "''L" •- t^-^^ : " , c3 r> ,1/ 

L .It 3,ljx, F3. ) , "/', f:.J, "/', fi. J ,13 X, fl. 0,2/ 

jt.-) M ,2 7 ,1U; /vC c b;./rt (..l;;vct C. '3 , ro-'l , r3 j , ro j , r3 7 ; wr t (. . o , F ^. ,.} ...' 

jrt J . 3 , .., 7; ./r t L,.2;-7tD !; ,2 7,10 

uit 3, 1-1 x, "'Ki'.3i 332 'Ot2 ", 2 X, "Pi-,? :-.if?", 3 X, "Vx/",11 X, "iin i 



'ic 



J 



JC C J. 3 

vru (,.■-; --^t":! '"■ ,^. 7 ,1.1; vr t u.Sr'vto t,27,l3;/'rt <:". 6; ,; t::' 5 ,2 7 ,1 J ; .vf. > j, 

L^ir T^l t.3 L3^;.vto (.,v5;i3Xt I;wtO G ,13 

: It 3,lx, ,;3. 1 ,3 X, f J. 2 ,3 X, f3. 1 ,3 X, f7. 1 , '. X , f 3. 3 ,j x, f 3 . I,'. >., t .d . 1 ,12 v, 

c It i , Ix, : ) . 1 , 3 X, t J . 3 , i 3 X, I J . i , J x, 1;3 . L , j x, fj . 1 

l; ,it 2,lx, rj. .i ,:; ><, fj . 2 ,! X, f3. 1 ,3 X, f/ . 1 ,3x , .P< . 3 , 3 x, i;3 . 1 

L n t ,S , i X, I > , 2 , Vx, c 3 . 1 , J -t, £ 7 . 1 , J X , f'J . 3 



-109- 



216: CM K':'>', t *. 2 

2.17: (: .\t 3 , Ix, 1) . 1 ,i >:, 1 1.'^ , 5 x, 13 . I , J x, T? . 1 

21J: e.it :,!>:, f';.l ,.? X, t J. 2 

21?: i'J Ell 7nto ■' T3l:t-?Ti^" 

220: ■' ', V : \ [1] *r2/; ', [ 11 +r2J; ': [ 11 ■*r2 

^2L: ':3r 1-2 to r35; (I-1)/2t23 



">?'? 




"ta^i:5":iE t 2-r4 = :i; o l*ex'iC od* p5/{ nl + rt) ) *ol ; jmp 2 
? }.*( ( ">2)-o3) /{ r?4f n3) ( "( P5* o5/( '14-02) ) *pl 
{ c2fp3) *( 1000/ oL) *p7- p3+p3f rGt 

"i3t\ot":if 0?r733";ito "store" 

236: l*I;if f:li7;r 50*r35; wrt C. 5 ,0 ,2 [ T ] + r30 , 1 [T ] ,'^ [I ,21 ; j 



? }.*( ( ">2)-o3) /{ r?4f n3) ( "( P5* o5/( '14-02) ) *pl 
{ c2fp3) *( 1000/ oL) *p7- p3+p3f rGt 
"i3t\ot":if 0?r733";ito "store" 

l*I;if f:li7;r 50*r35; wrt C. 5 ,0 ,2 [T] + r30 , 1 [T ] ,'^ [I ,21 ; j -i 
if El3lO;rfrt 6. o ,0 ,2 [ 1] + r3 J ? jmn 3 
vr t 6. 5,0 ,:[ 11 4t30, 1[T1 /"^f 1,2] 
if cl7 B;r 35+r53; r 50*c35 
£5r T=2 tor35;3ll 'ra tG'( 2 [ T ] ,2 [ I -1] ,: [ T ,11r^ f I .21 ,^ 

(T-1) /2*r23; i£ fig 7:imp 6 



■^ fT-1,2] ,z, ■;, V 



223 

22 3 

230: -11 'th 

231: itD ":> 
232 
233 
2 34 

237: - - - 

2 3 3: 

23 9: 

24 U: 

,n) 

241: (T-1) /2*r23; i£ fl3 7;jmp 

242: if f.Lg 10=0; j-io 2 

243: if >UI1=99 or fl^lJjsfg 10;JTir, 3 

24*: rfrt 6. 3 , Z-tc 33 ,0 [ l] ,5 [I ,1] ,W 

245: wrt 6 , r23 ,2 [ [ ] +r33 ,^ (I ] ,'^ (T ,21 , V," ,K [ I ] , L fl ] ;gto "rji" 

246: vrt 5. 4 , Z-IC30; wrt 6. 1 ,r23 ,2 [ I 1 + r30 /I ,K [ T ] , L [I ];gto "n 1 " 

247: wr t 5.3, Z-k:30 ,1 [I] ,E [I ,11 ,W 

248: wrt 6. 2 , r2 3 ,2 [ fl +r30 ,nI (H ,^ [I ,2] ,W ,4 

249: "ni":n3xt I 

250: -99*L[c35fl] 

251: gto "store" 

252: "rate ": (pi-p2i /2 + o2+p3; { o3- p5) /( p3-p2) *x>l 

253: (p4-p3) /( Dl-o3) +^6; ( pl-p2) /30*p9; ret 

254: if fig = 1; 7 to " ] u.no " 

255: "^tora":cl30 "D^TA OISC*"; sto 

255: "oj " :ar a C + D 

257: rS4 + I;r65*i; r 66+r; r 67+Mrf x-a 2 

258: 3tr (I ) *T $[ 1] ; T $[ 2] + I $f 11 

259: 3tr (S)+I$[3]; T$f 41* IS[ 3] 

250: str {T-1970) ^ T ?[ 5] ; T $1 6) * I $( 51 

261: ^Tar (int CVIOO) +61) * l$[ 61 

252: I$[l,ol*-l$ 

263: iso JI$, "correct =iat3 ^ ti-ne?";3to 

264: "a33a":3S3n f>'=;,l,0,< 

2.5 5: if X=l;jcD "0030" 

265: 3.irt 1 ,0 [ *] ,i: [ *1 , < [ *| , L [ *1 ,;^ $ ,3$ ,1 ,"; ,r , 1 ,1 , ^ , "en] " ;j ? ^ 2 

267: "3T3a":03en l'5,12;3t3 "asgn" 

263: Ejc3 Ujort T , V" , 5 , "/' , r , " ",q,"dat3 3t3r-J 13",I$ 

26 9: "32 ":if fl^ 5; 3 to "sf " 

270: if -l$!f"/23";3tp ydsp "Tress run" 

271: J3e3;l3p "set -plotter " ;3to ;if r22>7. 5 srd r23<i3r^t-:> 6,12;i-io 

272: :nt U,4/;.vct o 



27.i 


274 


t 


27 3 


2 73 


277 


Z7>; 


:'7D 


2o0 


^31 


:.;2 





J.17; 

■;:) J: 

2;h 

2''i: 

2 ■} 2 ; 

r :■ 3 ; 

2 "I ^r : 
2;i5: 
2."^'.; 
2 '17: 
•1 ■•>.•'. 

;i i-j -1 , 

3 J ^ : 
J ) 1 : 
T 

3.)2: 
J J 3 : 

T)7: 
}0 3: 

310: 
HI: 
212: 
213: 
Jil4: 
31 S: 
315: 
317: 
31-J: 

n): 

320: 
321: 
322: 
32 3: 
324: 
32 5: 
326: 
3^7: 
323: 
32:>; 



-110- 

: 3:1 -^,U.5,-.2,2.7;olt C,0,l;^lt ,. 1 ,2 

: lOC T-.l to 2.5 ;3y .l;->lt T,I;ilt -.2,T;olt .2,T;nlt U,T;n?>ct I;nl 
f ■• f 1 

: lit 1,0,2 

: ::>: T=l to 1 Jr ilt T ,0 

: if; T=ll or T = 2S;Ti.t T,.lr5lt I,api'.o I 

: -lit I,. 01 3; "Tit T ,-. L)13: -lit J ,0 

: ^2'<t T 

: JZiz 1.7S;:<1 i,:;olt lP,2.3,l;lol "^-s -:3?-:?9r ciiir" 

; -;i^ 1. Inj-Mt I'J ,2. -i 4 ,1; l^l r r 4 , "/", r .i 1 , "/' , r S5 

: Ut J i>,J. 3'M; I'll r f 7: "^^ 1 

: n t '1 , . 'j , It - -1 It -2.2,-. 3; 1 ^ 1 " . T" 

: -t: T-1 tci.j .;y . V; n t ),I,l;-'iIt -■^3,-.3;U)l T;n:'Xt J 

: 731.- 1.2:.;^lt l,,-,Jl,l;nt - . 'i^ ,. y^ 1 ,2; z-^lt - 2 . f ,- . 4 ; i ;l "30" 

: ^It ,m/l.") ) ),L; Ut -. ^ Sri )/i.)T ),2 r-^l t -2.5,-.4;lil "-7I " 

: 3 :; i ^ 1; c ^ 1 ' i 

; ci: T-2 tc 2J b/ ^'nt T/2,0,l;--^lt - ] . C ,- 1. 3; 1^1 T;n?xt T 
: M 1 , - 1 + 1: ^ : • , f i] >r 3 * r 1 ; r 3 5 + rn ; U -«• r 

:Tr I--2 t^ r 3 "); rin (r 2,MT ,11 ) -r^- ->in fr r,'MT ,21) -r>;">3:^t t 
= :>; "" ^2 torj3 

: '! ii (■.. r^: 1 ,r^) +r i; nT< f, fT 1 ,rt) *rl 
: i;' ("[r 1 +r3 J) /l,)n>- 2.'j1; [-l-..r21? iT-* 2 
1 : < ■: T 

c2-^rl3;i.C c 3<^2:3 or rl>^33)Mto "i;0-l^;j'' 

:>- T-1- to 21 :r/];->lt T,tM;7^1t -2,-1. 3;1M "^^^ + 2 r- ^;1'^.,t T 

" 1" 1-1^.1' : li .; + :- 

;">L- T=12 t.i M ■>/ ?;if '^='^2 0; 6C-*''" 

3lt T, 0,1— lit -2 ,-1. "!; L'nl Pf-^ + Cn + i-^n-xt I 

^Tt: T-v,j t-, n ^-v 2; Tit J, 0,1; colt -?,-1.2;j;t1 c ^ 3; r ^^ 3 + -^ *r n ; oo xt 

Tit J , ,1 , 1- : T 1 1 i. ,- 3: 1 T 1 "■,■; i n .1 , 3oee3 { r/ r-) " 

^It 13 ,..) ,1; : t1 t G,-3;1'tL "vin.l Hrection (3-;Tr3-:?T ) " 

>lt 2 .' ,.i ,1; ;ol t u,-3;1t1 " PDt 3 1 1 i t 1 Ti3nT . ( torj . C)" 

-3i3 1. 3,l,i,-0; Tit 0,. ^,1;--Tlt -.5,2,3:It1 "heirit t^-^-'-rl)" 

■: 1 i 7. . 7 'S , 1 , 1 , 1 

-; 5r r = 2 tor 21; IT K[ll = 'i rn3 L [ r 1 =!■; ] ct ■' r. f ] ■' 

Tit Un /2,(:[ r ) +r3 J) /I jj-); iolt u,0 

1"! Kt I 

"■".[1':t1c ll,0,l;Li; c3<=25 or i:i>=335;no +4 

■DC 1 = 2 tor21;ir: Ldl^'i r-n3 mi=j;Tt'T "v\]" 

Jit L fl 1/3M+12,(2[I 1 ^-r3 )) /n33 

i Tl t l: ,U; i?-<t T; JtT +7 

zd: 1-2 tor21;ie ^[11=0 ,^n.:i <(t]=0;3to "rf]" 

l f r, [T ] > = 1J,1; i -IT 3 

Tit L [T] /3J+lo ,(:[ [1 +r3 J) /I J]J 

i. T 1 1 U , 0; J t: T " .i 3 x t v i n i " 

Tit L[I]/3J + o ,(:[ [1 +r3 Ji /I ni; iolt 0,0 

"t3 vcn-i/inj " rnaxt I 

■•2[j':Tlt 25,J,l;if f 1 1 c = 0; r 35 -r53 

cor 1=1 toj:5J;if C [I ] +c3 J >25J3 ; gto "ott" 

if ELI ,11 =99 or 1 = 1; JTiT 3 

Tit (F[I ,l]-i:2) /2^25,(,(:[I] +::[I-1] )/2+r30) /I JOG 

iTlt 0,0; J[I ] +r3 4 

if L [I ,2] =y?; j no 3 

Tit (r;[I ,2) -r2) /24-23 ,(C[ r ] +r3j)/UJ.j 

iTlt 0,0; 2[t ] *r34 

M:Kt I 

" 3 t " : i £ 1 1 3 u; ] t T "2,\1" 



330: "^ij 



jju: 3 1 J •• : -111- 

3J2: int " wIt a O^ta ?" , v$; i f /^^jjip^fi 5 

33 3: if el} 5=0; ] CO "ort" 

334: or J=l tcr3-,;ie C[ 7 ] +r3 J > 1 ; J -1 *r5 3 ; j to "^val" 

33 o: V?v^l " :C*r5i>; O+rjJ 

33j: ^[i]3n{r.rj|)*(c(J]-':[J-l]} + r3J*r3J;next T 

339: •fJ].5ii ('. [T !)*(■.-:[ 7-11 -r3D+riJ*rn 

343: \(J|co3 {r,[] ])*{->:[J-i]-r3J) fr5:i*ri) 

^41: r b9/(,--r)j) *r3:i;c ( (./('i- ri3) *r5J; (r 5S *2+ r3:} ^2) \'^^v 

.42: pi^.M3tn{rS>/r,VJ)) >D;oct -^yif r>}>) an3 r.. j>j ,- i;o-^-. O; in-, 

jn: Lf : 59>0 ar.j r iO <J; ->+270 *0; j no ! 



1 L iT :)S-<0 



3nd r J JO; 270-1* r. jTin 2 



34 4; 

34 5: L •: : 5':.<3 fnl r 3 )> ]; i +<r;ri +0 

*i*.i "■^' '^" •'^^t J, j/ ,i J ;:,■'********* *****vv **v**,,t,, t,,, ,,,,,, ,,,,^^ ^^.^ 

347: e.Tt L,2J{ ,»• IVT].:. '^^T 3 IT ' , J x, f) . 1 , "i x ''^r^tcr-s " V 

34'^:^^l:nt W^Sx, "V-<yM x, H . 1 , 2 x, "V s ", 3 x, "n -fC i^IO i" , f m ,2 x, " le j .c . I . ^ 

S it' ;;U'l2-''^^'' fl-3r;,vrt C.2;.vrt C.3,.,.;.rt G 

3!'.3: ;t3 
* 2 74 12 



-112- 



Appendix 3 
Hourly Average Solar Radiation Data 
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Appendix 'f 

Potential Temperature Profile Plots 

May to September, 1977 to 1979 
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Appendix 5 
Solar Flux Analyses 

Section A5(a); All days 
Section A5(b): Convective Days Only 
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Section A5(a) 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 

Month: JANUARY 1979 

Solar Flux (Watts/Sq .M. ) 

Standard # of 

Hour Maximum Minimum Average Deviation Data Points 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 

Month: FEBRUARY 1979 

Solar Flux (Watts/Sq.M.) 

Standard # of 
Hour Maximum Minimum Average Deviation Data Points 

28 

Q 23 

23 

28 

5 9 25 

16 77 53 25 

57 190 97 25 

49 297 133 25 

66 373 153 25 

49 402 166 25 

57 400 154 . 25 

74 349 130 27 

41 275 121 28 

25 153 ?4 28 

8 3 3 21 23 

Q Q 28 

a i 28 

a Q 23 

28 

Q 23 

Q 23 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 

Month: MARCH 1979 

Solar Flux (Watts/Sq .M.) 

Standard # of 
Hour Maximum Minimum Average Deviation Data Points 

31 

31 

31 

2 11 30 

54 65 30 

140 113 30 

33 239 158 30 

57 347 20 3 30 

57 449 230 30 

41 518 241 30 

66 511 2 31 30 

57 437 213 30 

41 32 5 175 30 

8 205 124 30 

80 6 9 30 

8 24 30 

31 
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31( 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 



Month: APRIL 1979 



solar Flux (Watts/Sq .M. ) 

Standard # of 
Hour Maximuiti Minimum Average Deviation Data Points 
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t 30 

n 30 

7 3 4 7 20 

M 184 111 20 

ii 319 183 20 

M 4 36 2 31 20 

10 532 26 20 
t$ 59 277 20 
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11 58 2 32 20 
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» 268 158 19 
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i 28 28 18 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 

Month: SEPTEMBER 1978 

Solar Flux (Watts/Sq .M. ) 

Standard # of 

Hour Maximum Minimum Average Deviation Data Points 

30 

30 

30 

30 

25 21 30 

106 85 30 

8 222 155 30 

25 294 200 30 

33 373 240 30 

33 415 275 30 

16 431 284 30 

25 403 266 30 

25 350 212 30 

16 265 184 30 

8 166 123 30 

6 5 53 30 

9 14 30 

30 

30 

30 

30i 



1- 4 





4- 5 





5- 6 





6- 7 





7- 8 


74 


8- 9 


246 


9-10 


451 


10-11 
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11-12 


721 


12-13 
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13-14 


901 


14-15 


836 


15-16 


713 


16-17 


565 


17-18 


360 


18-19 


205 


19-20 


49 


20-21 





21-22 





22-23 





23-24 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 



Month: OCTOBER 



1978 



Hour 



Solar Flux (Watts/Sq.M,) 

Standard # of 
Maximum Minimum Average Deviation Data 



Points 



1- 4 





4- 5 
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5- 6 
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6- 7 
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7- 8 


82 


8- 9 


26 2 


9-10 


393 


10-11 


475 


11-12 


573 
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606 
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14 7 
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294 

324 
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306 

244 

159 

80 

12 
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18 

57 

104 

14 7 

158 

187 

199 

194 

172 

123 

76 

24 













31 
31 
31 
31 
31 
31 
31 
31 
31 
31 
30 
31 
31 
31 
31 
31 
31 
31 
31 
31 
31 
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MINISTRY OF THE ENVIRONMENT 
SUDBURY ENVIRONMENTAL STUDY 

SOLAR FLUX ANALYSIS: 



Month: NOVEMBER 1978 

solar Flux (Watts/Sq .M. ) 

Standard # of 

Hour Maximum Minimum Average Deviation Data Points 

30 

30 

30 

30 

9 10 30 

64 51 30 

14 9 92 30 

16 222 125 30 

16 264 152 30 

16 26 5 154 30 

8 229 134 30 

165 101 30 

84 66 30 

22 29 30 

1 2 30 

30 

30 

30 

30 

30 

30 
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